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ABSTRACT 

The nucleus of M87 displays a LINER spectrum at optical wavelengths, with a nu- 
clear disk of nebulosity that is resolved by the Hubble Space Telescope. We present 
new results from optical and ultraviolet spectra of the central ~ 40 pc as measured by 
HST. In contrast with previous results for the extended disk, the emission-line spec- 
trum of the central region is best reproduced by a multi-component photoionization 
scenario, rather than shock heating. The nebular properties as well as energetic con- 
siderations suggest a transition on scales of several tens of parsecs, from predominantly 
photoionization by a central accretion source, to shock ionization within the disk. If 
this source is at all representative, it suggests that many LINERs may be composite in 
terms of the energetic processes that give rise to the emission spectrum. We also report 
measurements of resonance-line absorption for the nucleus. The absorption spectrum, 
like the emission lines, is characterized by low ionization. The absorption line measure- 
ments coupled with independent constraints suggest a total hydrogen column density 
of 10 19 — 10 20 cm -2 , outflowing from the galaxy center with a velocity of ~ 126 km s . 
The kinematic signature of an outflow, along with evidence that the absorber covers the 
power-law continuum source but not the emission-line clouds, suggests that the absorb- 
ing matter is related to accretion phenomena in the nucleus. The presence of such an 
outflow resembles similar behavior in luminous AGNs, although the low ionization that 
characterizes LINERs is probably indicative of a different mode of accretion in these 
sources. 

Subject headings: galaxies: individual (M87) — galaxies: nuclei — galaxies: active 
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1. Introduction 

A large fraction of nearby galaxies harbor low-ionization nuclear emission-line regions (LIN- 
ERs; e.g., Heckman 1980; Ho, Filippenko, & Sargent 1997). A variety of mechanisms have been 
suggested for the underlying power source in these objects, including shocks, black hole accretion, 
and hot stars. LINER case studies have revealed candidate prototypes for each of these phenomena, 
strengthening speculation that LINERs are a heterogeneous class (e.g., Filippenko 1996). 

The LINER galaxy M87 is a particularly interesting example. Ground-based optical spectra 
of its nucleus show narrow-line emission with a complex velocity profile, but no clear signature of 
the broad permitted lines that are the hallmark of quasar activity. Observations with the Hubble 
Space Telescope (HST) reveal a spatially resolved disk of emission at the nucleus, and analysis of 
emission from the disk provides strong evidence of shock excitation (Dopita et al. 1997, hereafter 
D97). Considered alone, these properties give little indication that the system is related to luminous 
AGNs; yet this object provides one of the best cases for accretion-powered activity in the nucleus of 
a galaxy. The evidence takes the form of its famous synchrotron jet (e.g., Boksenberg et al. 1992; 
Bicknell & Begelman 1996) and a velocity field in the nuclear gas disk requiring the presence of a 
central dark object with a mass of 3.2 x 10 9 M (Ford et al. 1994; Harms et al. 1994; Macchetto 
et al. 1997). 

One possible interpretation of the existing results for M87 is that LINER emission is caused 
primarily by shocks, even when coexistent with an accretion power source that may generate sub- 
stantial ionizing radiation. This scenario contrasts markedly with the behavior of Seyfert galaxies, 
which have narrow- line emission powered primarily by photoionization (e.g., Laor 1998). Moreover, 
the evidence for shocks in the M87 disk comes from ultraviolet emission-line ratios, but the same 
diagnostics observed in a handful of other LINERs do not provide similar indications of mechanical 
heating (e.g., Barth et al. 1996, 1997; Maoz et al. 1998). The latter finding may imply that the 
LINER in M87 is an anomaly; however, another possibility is that the uncharacteristic result for 
this source reflects the unusual circumstance that its central nebulosity is resolved, and the previous 
HST study considered emission from only a part of the disk that did not include the nucleus per 
se. In order to understand fully the nature of the nebulosity at the center of M87, and to assess the 
larger implications for LINERs, measurement of the UV/optical spectrum of the nucleus is clearly 
desirable. In this paper we present the results of spectroscopic observations carried out with HST 
for this purpose. 



2. Observation and Data Analysis 

Long-slit spectra of the nucleus of M87 were obtained with the Space Telescope Imaging 
Spectrograph (STIS) on 1999 February 1 UT, using the 52" x Of! 5 aperture and the gratings listed 
in Table 1. The slit was oriented at position angles 216° and 240° for the optical (CCD) and 
UV (MAMA) data, respectively. (The difference in position angle resulted unintentionally from a 
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roll-angle adjustment during the observation sequence.) The resulting spectra provide continuous 
wavelength coverage of the nucleus in the range ~ 1150 — 10270 A. The relatively large slit width 
leads to rather low spectral resolution (see Table 1) but was selected in order to maximize the 
measured signal in the UV bandpass, while maintaining a uniform aperture dimension across the 
entire spectrum. The full-width at half maximum (FWHM) spatial resolution is ~ / /05 for the 
MAMAs, and ~ / . / 12 for the CCDs. The data were calibrated "On-The-Fly" upon retrieval from 
the data archive at the Space Telescope Science Institute. All data reduction and analysis was 
performed with IRAF 8 . With this slit width, geocoronal Lya is blended with the redshifted Lya 
from M87. We measured the background geocoronal emission in the G140L spectrum near the ends 
of the slit, and subtracted the result from each row of the two-dimensional spectrum, making sure 
that we also propagate errors in the correct way. 

Spectra of the nucleus were extracted for an aperture size of ~ 0?25 along the slit (i.e., 10 
MAMA pixels and 5 CCD pixels, corresponding to approximately twice the spatial FWHM for 
the CCD). The result is shown in Figure 1. In order to aid in removing any underlying galaxy 
starlight, we generated a stellar continuum template by averaging the spectra of off-nucleus regions 
that were emission-free. The resulting aperture was 10" (5" for UV spectra) wide and centered at 
7'.'5 (3'.'75 for UV spectra) southwest of the nucleus. Emission from the nebular disk at these scales 
is negligible within our aperture, in comparison with that measured from the nucleus. HST images 
show a rapid fall-off of nebular surface brightness with radius (e.g., D97). In contrast with the FOS 
study reported by D97, which measured the spectrum tf.'G from the nucleus with a l'.'O aperture, 
the data reported here employ a narrower slit, and consequently do not exhibit a signal-to-noise 
ratio sufficient to measure nebular features outside the central aperture defined above. 

We measured the fluxes of the emission features while simultaneously fitting the continuum, 
via x 2 minimization, using SPECFIT as implemented in IRAF (Kriss 1994). We fitted Gaussians 
to the emission lines, and represented the continuum as the sum of a featureless power law and 
the galaxy template. The results of the measurements are given in Table 2. In all cases, the fit 
results were consistent with the continuum being dominated by the power-law component. Errors 
were evaluated by SPECFIT from the accompanying error extensions to the spectra. The FWHM 
velocity for H/3 is ~ 2100 km s -1 , after correcting for the instrumental profile for a point source, 
and ~ 900 km s _1 if the emission is extended uniformly across the slit; the latter value can thus be 
taken as a lower limit. Due to the significant intrinsic velocity width and the low spectral resolution, 
Ha was heavily blended with [N II] AA6548, 6583, and the [S II] AA6717, 6731 lines were also heavily 
blended with each other; no attempt was made at deblending these composite features, which we 
instead fit by single Gaussians. The width and redshift of the [O I] A6364 template were linked to 
those of [O I] A6300 and, as dictated by atomic physics (Osterbrock 1989), [O I] AA6300/6364 = 3.0. 
The same procedure was followed for [O III] AA4959, 5007. 



8 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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A joint analysis of the optical and UV spectra is potentially sensitive to reddening. The fact 
that Ha was heavily blended with the [N II] lines makes it impossible to get a direct measure of 
reddening using the Ha/H/3 ratio. Moreover, H5 is blended with [S II] AA4069, 4076, and H7 leads 
to unphysical extinction values when combined with H/3; this may be due to an overestimate of the 
H7 flux caused by blended emission in [O III] A4363. Our spectral range includes the recombination 
lines He II A4686 and He II A1640. For Case B recombination, the ratio He II AA1640/4686 » 7 
(Hummer & Storey 1987). Our spectra set a 3a upper limit on He II A4686 of 20% of H/3. This limit 
is calculated according to the formula 3\/27rcr c (14A), where a c is the root-mean-square uncertainty 
per Angstrom in the local continuum, and 14 A is the measured Gaussian dispersion of the H/3 
profile. (This prescription corresponds to a Gaussian profile with amplitude 3cr c .) The result is 
(He II AA1640/4686) o i> ser?)e d > 2.0. Employing the extinction curve of Cardelli, Clayton, & Mathis 
(1989), we find that Ay < 2.4 mag. The Galactic foreground extinction toward M87 is expected to 
be much less than this value, and independent estimates of the total absorption (§3) suggest that 
the Galactic component dominates the total extinction, with Ay « 0.15 mag. 

Additional spectra of the nucleus of M87 spanning 1140 - 6820 A were obtained on 1996 July 
2 UT with the FOS, using the 0'. / 21 aperture (Table 1). The calibrated data were obtained from 
the public domain archive. We measured the emission lines as for the STIS spectra, assuming a 
power-law form for the continuum. Uncertainties in the line fluxes are generally larger for the FOS 
observations, reflecting a somewhat lower signal-to-noise ratio. Measurements from the two spectra 
are in general agreement although significant differences are present in Lya, which may be affected 
by uncertain geocoronal subtraction for the FOS spectrum, and [Ne III] A3869. 

In addition to emission lines, there is a wealth of UV absorption features seen in the spec- 
tra. The velocity shifts of some lines are consistent with absorption within the Galaxy, while 
others appear associated with M87. We used spectra with the highest spectral resolution for our 
absorption-line measurements. Table 1 shows that all the FOS gratings have higher spectral reso- 
lution than those employed with STIS, except for FOS G160L, which has a lower resolution than 
STIS G140L. The absorption lines in the FOS spectra segregate into two systems distinguished by 
their line widths; those with FWHM k, 200 km s _1 are of Galactic origin (and are unresolved), 
while those with FWHM ps 400 km s _1 are intrinsic to M87, where the FWHM values quoted 
here have not been corrected for the instrument profile of FWHM= 230 km s _1 . The absorption 
lines in the STIS G140L have FWHM ~ 400 km s _1 , uncorrected for the line-spread function of 
FWHM = 190 km s" 1 for a point source; the velocity widths of many features in this spectrum are 
potentially affected by blends. Identification of features was based on wavelength correspondence 
with line lists from Verner, Barthel, & Tytler (1994), and consistency tests with the measured line 
strengths or upper limits for resonance features from the same ions. 

The STIS and FOS observations are separated by ~ 2.6 years. The equivalent widths of the 
absorption lines agree to within the measurement errors. The continuum level in the FOS spectrum 
is higher by 20%, with greater contrast at shorter wavelength. Tsvetanov et al. (1998) discussed 
variability in the optical continuum during 1994 and 1995, and found differences of a factor of ~ 2 
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on a timescale of 2.5 months, so variability is clearly a concern when comparing measurements from 
different epochs. 

3. Results and Interpretation 

The physical interpretation of nebular emission in LINERs is complicated by the fact that 
several different energy sources may give rise to such emission. Theoretical predictions of optical 
lines are degenerate in the sense that the same line ratios can be explained equally well by more 
than one physical scenario, including in particular shocks and photoionization (e.g., D97; Allen, 
Dopita, & Tsvetanov 1998). Allen et al. (1998) have suggested that the inclusion of ultraviolet 
emission diagnostics provides a basis for distinguishing between shocks and photoionization, and 
UV emission lines figure prominently in the analysis by D97 implicating shocks in the extended 
nuclear disk in M87. In the present study we extend this analysis to the interpretation of the 
emission of nuclear emission as measured in the central ~ Of! 5, corresponding to ~ 39 pc (for an 
assumed distance of 16.1 Mpc; Tonry et al. 2001). 

3.1. Line Ratio Diagrams and the Excitation Mechanism 

Comparisons of our measurements with theoretical predictions are presented in two-dimensional 
line ratio diagrams in Figure 2, which employs optical and UV lines, and Figure 3, which uses ex- 
clusively UV lines. These diagrams are selected from those shown by D97 and Allen et al. (1998) 
to hold promise for distinguishing ionization processes. In addition to our STIS results, line ratios 
from the FOS spectra of the nucleus (our measurements) and of the ionized disk (D97) are also 
shown. The diagrams involve intensity ratios of lines close in wavelength, so that they are rela- 
tively insensitive to reddening; the plotted data points have consequently not been corrected for 
reddening. 

Predictions for several excitation mechanisms are overplotted in Figures 2 and 3. Predicted 
ratios for shock excitation are taken from Dopita & Sutherland (1995); for their models, the pa- 
rameters that produce the range in line ratios are the shock velocity and the magnetic parameter 
B/n 1 ! 2 , where B is the magnetic field strength and n is the hydrogen number density. Predic- 
tions for photoionized plasma assume power-law continua with a spectral index a = (—1.0, —1.4) 
(assuming f v oc v a ), density nn = (100, 1000) cm" 3 , and ionization parameter U = 10~ 4 — 10 _1 , 
where U is the ratio of ionizing photon and hydrogen densities at the face of the irradiated cloud 
(Binette, Wilson, & Storchi-Bergmann 1996). 

A possible limitation of these photoionization models is their assumption of a single population 
of ionization-bounded clouds. Binette et al. (1996) have argued that the spectra of AGN narrow- 
line regions can be better reproduced by photoionization scenarios with multi-component cloud 
systems; in particular, they emphasize the possible importance of a mix of ionization- and matter- 
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bounded clouds. We have consequently also shown for comparison in Figures 2 and 3 the predictions 
of a two-zone photoionization model, in which we have reproduced the Binette et al. results using 
the photoionization code CLOUDY, version 90.05 (Ferland et al. 1998). In this scenario, matter- 
bounded (MB) clouds filter the ionizing continuum from the AGN before it irradiates the ionization 
bounded (IB) clouds. The MB clouds are irradiated by a power law with a = —1.3, and arc 
characterized by Umb = 0.04 and hub = 50 cm -3 , while the IB clouds have U = 5.2 x 10~ 4 
and tijb = 2300 cm -3 . The relative contribution of emission from the two cloud populations is 
parametrized by Am/i, which is the ratio of solid angles subtended by MB and IB clouds. Formally, 
A m ij is always larger than one, but in calculating the resulting spectrum, it can be effectively less 
than one if the MB clouds are hidden along our line of sight but seen by the IB gas. 

The measured ratios reported by D97 at an off-nuclear position, and for the nucleus as given 
here, are different, and Figures 2 and 3 suggest different interpretations for the two locations. A 
similar inference was noted by Sankrit, Sembach, & Canizares (1999) based on FOS measurements 
of the C IV/Lya ratio alone. As discussed by D97 and illustrated also here, the line ratios for the 
off-nucleus measurement of M87 are consistent in all of these diagrams with a shock interpetation; 
the photoionization predictions, particularly the single cloud models, do not agree with the data. 
The line ratios measured for the M87 nucleus are, however, less clear-cut in their interpretation. 
In Figure 2 (optical and UV lines), the nucleus ratios generally fall within the shock model loci 
and are inconsistent with both single-cloud and A M n photoionization predictions. However, in 
Figure 3 (UV lines) the nucleus measurements fall outside the predictions for shocks. Single-cloud 
photoionization models are likewise inconsistent with these data; but the A^n sequence intersects 
the measured ratios in all of the plots in Figure 3, and in most cases the data are in fact consistent 
with predictions for a single value of A^n ~ 0.03. This result provides a strong suggestion that 
photoionization is, in fact, important in the nucleus. In contrast, comparison of the observed ratios 
with the detailed shock models of Dopita & Sutherland (1995) shows that the optical line ratios 
in Figure 2 are not reproduced by a single combination of shock parameters. The A M /j model 
remains problematic, however, for explaining the optical line ratios in Figure 2 as well as generally 
underpredicting the strengths of the UV lines relative to H/3 for the nucleus. 

Is it possible to reconcile a single excitation scenario with the nucleus line ratios plotted in both 
Figures 2 and 3? It is interesting to note that the division between optical and UV lines separates 
strong transitions not only by excitation energy, but also by the critical densities for the collisionally 
excited lines: lower n cr it for the optical lines, which are forbidden transitions, and higher n cr n for 
the UV lines, which are intercombination or resonance features. As a result, increasing the plasma 
density will affect the predicted line strengths in different ways; the forbidden lines will become 
suppressed, while the lines with high n cr n will strengthen as they take over a larger share of the 
nebular cooling. For the M87 nucleus, increasing the nebular density might thus be expected to 
improve the success of photoionization models by decreasing the [O IIl]/H/3 ratio and boosting the 
strengths of the UV lines relative to H/3. Nebular components spanning a wide range of density 
are known to exist within AGNs, as revealed by line width-n cr jt correlations (e.g., Pelat, Alloin, 
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& Fosbury 1981; Filippenko Sz Halpern 1984; Filippenko 1985), and dense components on small 
spatial scales are explicitly revealed in nuclei through HST studies (e.g., Barth et al. 2001). 

To investigate these effects, we experimented with A M /j sequences with increased densities 
and found that many characteristics of the nuclear spectrum are reproduced if hmb = 10 6 cm" 3 
and njB = 10 6 ' 3 cm" 3 . The calculation results are again shown for comparison in Figures 2 and 3. 
The observed line ratios fall consistently near the predictions for the high-n A M /j 0.002 — 0.004, 
except in Figure 2a. In the latter case, the predicted [O II] A3727 strength is very weak, due to the 
low critical density of this line. A likely remedy to this problem is the inclusion of an additional low- 
density IB (low-n IB) component that emits efficiently in this transition while minimally perturbing 
the other line ratios. 

We modeled this third component as the low-density IB of Binette et al. (1996), except that 
it was irradiated by an ionizing continuum, the same power-law continuum discussed above, which 
has not been filtered through the MB clouds. To achieve the final combination of high-n A^n + 
low-n IB, we assumed that our observation aperture is filled with high-n A M /j = 0.004 nebulae and 
low-n IB clouds. Relevant line ratios for each component, normalized by H/3, are listed in Table 
3. The composite model represents a weighted sum of the three components, and a good overall 
match with the nucleus observations was obtained with relative H/3 contributions from (high-n 
MB): (high-n IB): (low-n IB) = 0.02 : 0.68 : 0.3 (note that the high-n MB : high-n IB ratio is 
already fixed by A m ij). The composite predicted ratios and observed values are listed in the last 
two columns of Table 3. 

While the 3-component simulation should not be taken too literally as a physical representation 
of the nebular structure in the M87 nucleus, the comparison presented here nonetheless provides 
an important demonstration that a multi-component photoionized plasma may account for the 
emission-line properties of this source. The choice of nebular parameters in this picture is not 
entirely arbitrary. Gas with densities ranging from ~ 100 cm" 3 to ~ 10 6 cm -3 exists in the cores 
of AGNs, and more general arguments exist for the importance of composite cloud populations in 
these environments. In particular, Baldwin et al. (1995) and Ferguson et al. (1997) have discussed 
the phenomenon of locally optimally emitting clouds (LOCs), in which emission in a particular line 
emerges predominantly from the subset of clouds with appropriately favorable nebular conditions 
for that transition. In the spirit of the LOC models, we conclude that it appears likely that a 
combination of high-density nebulae emit most of the lines observed in the nucleus of M87, while 
lower density clouds lead to the emission from ions with low critical densities, such as [O II] A3727 
and [S II] AA6717, 6731. Analyses of other LINERs have similarly suggested that multiple cloud 
populations may be required to account for the observed emission- line properties (e.g., Pequignot 
1984; Barth et al. 2001). 

An interesting point which our high-n Aj^/j calculations raise is whether the UV line-ratio 
diagrams are in fact robust in distinguishing between shock heating and photoionization, if density 
effects are taken into consideration. Figures 2 and 3 show that the line ratios from the ionized disk 
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of M87 could be well described by the high-n A M/I = 0.006, except for [O II] A3727. We again 
attempted to construct a composite model with different combinations of the high-n A^/j and 
low-n IB components. We found that the same fractional contributions of the components used 
above (with H/3 proportions 0.7 : 0.3) give reasonably satisfactory results, although low-n cr it lines 
such as [O II] and [S II] AA6717, 6731 remain underpredicted by a factor of ~ 2. We conclude that 
the extended disk of M87 remains a good candidate for shock ionization, but emphasize that the 
UV line diagnostics advocated by D97 and Allen et al. (1998) may encounter limitations when 
applied to composite nebular systems that include high-density components. 

3.2. Energetics 

A consistency test for photoionization interpretations of the nuclear nebula in M87 is whether 
the central source provides sufficient ionizing photons to power the line emission. To address this 
point, we employed recent measurements of the nuclear continuum as measured by XMM-Newton 
in 2000 (Bohringer et al. 2001). These data have the advantage of simultaneously sampling the 
UV and X-ray bandpasses; while the central source is known to be variable, we employ these 
measurements as a representative snapshot of the broad-band continuum. The narrow lines in this 
source are neither known nor expected to be variable. In the UV bandpass, the XMM-Newton data 
provide a flux density for the M87 nucleus of (1.90 ± 0.04) x 10" 27 erg s" 1 cm" 2 Hz" 1 at 2120 A, 
and an integrated X-ray flux of 1.5 x 10 -12 erg s" 1 cm" 2 for 2-10 keV. 

We can extrapolate the XMM-Newton UV flux by assuming that the spectral shape is the 
same as during our HST observation. We first attempted to derive an analytic fit to the continuum 
we measured for the nucleus, using SPECFIT as before. The continuum was represented by a 
broken power law subject to reddening, combined with a stellar component represented by the 
template described in §2. The fit indicates that the galaxy starlight contribution is negligible, and 
the continuum is well-represented by only a power law with a break at 4501 ± 6 A, subject to 
reddening given by Ay = 0.15 ± 0.01 mag. This amount of reddening is reasonably consistent 
with the value of 0.11 ± 0.02 mag predicted from the Galactic H I column density toward M87 
[N H i = (2.1 ±0.3) x 10 20 cm" 2 ; Sankrit et al. 1999] and a standard N HI /A V ratio (Bohlin, Savage, 
& Drake 1978). The spectral indices of the power law are a = -1.35 ± 0.02 and a = -1.78 ± 0.02 
for the red and blue portions, respectively. The spectrum of the nucleus along with this fit are 
shown in Figure 4. A similar fit to the FOS spectra was reported by Tsvetanov et al. (1999b). 
Our 1999 STIS flux at 2120 A is (7.4 ± 0.6) x 10" 28 erg s" 1 cm" 2 Hz" 1 , a factor of 2.6 below 
that obtained in 2000 by XMM-Newton, which provides one indicator of the level of variability 
between these two measurements. Scaling our continuum observation to bring it into agreement 
with the XMM-Newton data, and removing reddening corresponding to Ay = 0.15 mag, implies a 
flux density at 2500 A of 3.6 x 10" 27 erg s" 1 cm" 2 Hz" 1 . 

At higher energies, the continuum is constrained by the XMM-Newton X-ray measurements. 
Bohringer et al. (2001) reported that the 0.6 — 10 keV spectrum was successfully described assuming 
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Galactic absorption only, and a power law with best-fitting a = —1.2. Using the same spectral 
index with the 2 — 10 keV flux noted above implies that the 2 keV flux density is 2.2 x 10 -30 erg 
s _1 cm -2 Hz -1 . The two-point spectral index connecting 2500 A and 2 keV is then a ox ~ —1.2. 
While this value is the same as that measured in X-rays, it is shallower than the measured UV 
index of —1.78, implying that the spectrum must show some spectral curvature across the ionizing 
ultraviolet region. If we approximate the continuum with a single power law with a f» —1-2, 
normalized to match the dereddened 2500 A flux, the hydrogen-ionizing photon flux at the Earth 
in the absence of absorption would be 0.12 photons s -1 cm -2 . After correction for Ay = 0.15 
mag, the H/3 photon flux in our STIS nucleus aperture is 1.90 x 10 -3 photons s^ 1 cm~ 2 , implying 
a photon ratio Qion/Qu/3 = 63. For comparison, Case B recombination at a temperature of 10 4 K 
and density of 10 3 cm~ 3 predicts Qion/Qup = 8-55 (Hummer & Storey 1987). It thus seems likely 
that the central source provides enough ionizing photons to power the emission nebula within the 
central ~ 1", if the covering factor for the nebular gas is ~ 10 — 20% . 

The nebular disk extends beyond the aperture employed here for the nucleus, however. The 
surface brightness distribution for the nebular disk has been studied by Ford et al. (1994), D97, 
and Tsvetanov et al. (1999a) using HST narrow-band imaging of Ha + [N II]. Within 1" of the 
nucleus, the resulting integrated Ha + [N II] flux is (2.0 ± 0.7) x 10~ 13 erg cm~ 2 s _1 (Ford et al. 
1994). If we integrate out from this radius the analytic fit provided by Tsvetanov et al. (1999a) 
for the surface brightness distribution, the total extrapolated Ha + [N II] flux for the disk is then 
2.6 x 10~ 13 erg s _1 cm~ 2 . This flux exceeds our nuclear measurement by a factor of ~ 3.3; if we scale 
our HP flux by a corresponding amount, the ratio Qi on /Qu/3 ~ 19, which still allows consistency 
with photoionization if the covering factor is at least ~ 40 — 50%. 

The calculation described above for the ionizing photon flux is uncertain in several ways. First, 
it requires an uncertain extrapolation across the unobservable ionizing UV bandpass that contains 
most of Qion- In addition, the continuum source is variable. It is therefore uncertain whether 
our choice of flux level adopted for this calculation is close to an average value. The calculation 
above involves a further assumption that the continuum we measure is emitted isotropically at 
the source. If a portion of this continuum is in fact beamed in our direction, our estimate of the 
photon production available for photoionizing the nebula would be overestimated. Some amount 
of beaming would not be surprising in light of other similarities between the M87 nucleus and BL 
Lac objects (see Tsvetanov et al. 1998 for discussion). In any case, it appears likely that the 
energetics change from being driven primarily by the nuclear radiation field, to other processes at 
larger radii; on large scales (< 10"), the total integrated line flux grows to 1.2 x 10~ 12 erg s _1 
cm -2 , implying as above Qion/Qnp ~ 4, which is likely to be problematic. Taken at face value, the 
models for photoionization of the nucleus (§3.1) and shock heating of the disk (D97) suggest that 
this transition occurs within only a few tens of parsecs from the center. 
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3.3. Absorption Lines 

The absorption lines detected in our spectra provide another diagnostic of material that may 
be in close proximity to the AGN. The data with the best combination of (high) spectral resolution 
and signal-to-noise ratio, from STIS G140L and FOS G270H, can be used to estimate some of the 
properties of the absorbing matter. Absorption lines identified in these spectra and measurements of 
their equivalent widths (EWs) are presented in Table 4, where we have also included measurements 
of the Ca II K + H and Na I D lines from the FOS G400H and G570H spectra, respectively. The 
lines were measured using SPECFIT assuming Gaussian profiles. We show in Figure 5 close-up 
views of the spectra and fits around the absorption lines. Tsvetanov et al. (1999b) have reported 
the detection of absorption lines in the same FOS spectra that we are using. 

The velocity of the absorption lines relative to that of the emission lines in M87 is of importance 
for determining the nature and physical state of the absorber. The measured wavelengths listed in 
Table 4 reflect the wavelength scale resulting from the pipeline reduction of the spectra, which may 
be subject to significant zero-point shifts. The average velocity of the absorption lines measured 
with the FOS G270H, not including lines blended with Galactic features, is 1091± 11 km s _1 ; for the 
STIS G140L spectrum, the mean velocity is 1066 ± 36 km s _1 . These values are consistent within 
their uncertainties, and for each grating considered individually the velocities of the individual lines 
are likewise consistent with each other. We also fit emission lines in the same spectra. For FOS 
G270H, C II] A2326 appears at 2334.9 ± 0.5 A, corresponding to a velocity of 1073 ± 66 km s _1 ; 
for STIS G140L, C IV A1549 appears at 1555.3 ± 0.6 A and He II A1640 appears at 1647.0 ± 0.6 
A, yielding velocities of 1220 ± 120 km s _1 (assuming an optically thin C IV doublet ratio) and 
1200 ± 110 km s" 1 , respectively. These results indicate that the absorber and the emitter are 
at approximately the same redshift. The G140L measurement yields a wavelength for Lya of 
1221.45 ± 0.06 A, corresponding to 1426 ± 15 km s -1 , which is significantly larger than the results 
for the other lines; however, this finding is consistent with expectations that the Lya emission may 
be modified significantly (30% - 40% flux reduction) due to absorption by Galactic H I, resulting 
in a redward shift of the transmitted line centroid (Sankrit et al. 1999). 

We can improve estimates of absolute velocity from these spectra by associating the Galactic 
absorption features with H I 21 cm absorption measured toward M87. Results from Davis & 
Cummings (1975) indicate that radio absorption arising in the Milky Way along this sightline is 
spread over a 50 km s _1 range, with the major component at ~ —8 km s _1 with respect to the Local 
Standard of Rest; this velocity corresponds to a heliocentric value of ~ —12 km s -1 . The average 
redshift of Galactic lines in the FOS G270H spectrum is nominally —72 ± 16 km s -1 , implying 
that a correction of +60 ± 16 km s _1 should be added to velocities derived from this spectrum in 
order to obtain heliocentric values. The heliocentric velocity of the M87 absorber is then 1151 =L 19 
km s _1 , based on the G270H measurement. This result means that the absorber is blueshifted by 
~ 126 ± 19 km s _1 with respect to M87, which is at 1277 ± 2 km s _1 as determined from stellar 
absorption lines (van der Marel 1994). The data for C II] in the G270H spectrum likewise implies 
that the emission-line plasma is also blueshifted with respect to the host galaxy, by 144±68 km s _1 . 
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We carried out a curve-of-growth analysis to calculate the column densities of the absorbing 
ions. We followed the procedure outlined in Spitzer (1978) using our observed equivalent widths. 
One major uncertainty in this derivation is the detailed shape of the absorption features, since the 
lines that we detect are not entirely resolved; some substructure within the lines is suggested from 
previous observations of the Na I D lines (Carter, Johnstone, & Fabian 1997). We consequently 
chose a Doppler parameter b = FWHM/1.665 = 200 km s" 1 for consistency with the intrinsic 
FWHM of the lines as measured in the FOS and STIS spectra. We used the oscillator strengths 
from Verner et al. (1994). The resulting values of the column densities Ni on are listed in Tables 4 
and 5. Given the low spectral resolution and possibility of unresolved substructure in the lines, the 
resulting N,i on estimates should be treated as lower limits. This point is reinforced by the tendency 
for some ions to display lower calculated column densities for lines with larger oscillator strengths, 
consistent with the presence of saturated substructure within these lines. 

The results for the individual ions can be used to place constraints on the total column density 
of the absorber. Each Ni on gives a lower limit to the corresponding N e i ement , which in turn provides 
a lower limit to the total hydrogen column density Nh, assuming solar abundances (Anders & 
Grevesse 1989). The inferred Nh from the various ions, whose absorption lines are not contaminated 
by absorption from other absorption features, indicates that Nh > 10 18 cm" 2 , except from the 
column densities of Ca II (N H > 10 19 cm" 2 ) and Mn II (Nh > 10 20 cm" 2 ). 

The total Nh in the M87 absorber can be constrained by other means. Dwarakanath, van 
Gorkom, & Owen (1994) used H I 21 cm observations to place an upper limit on the absorbing H I 
column density of 5 x 10 19 cm -2 . (This value is sensitive to the spin temperature, assumed by these 
authors to be 100 K.) Since the absorption lines detected in the HST spectra arise only from very 
low ionization species (requiring ionization energies < 13.6 eV to produce), the hydrogen in the 
absorber is expected to be mostly neutral, and this value can thus be taken as a useful limit on the 
total Nh- As noted in §3.2 for the nucleus, the reddening of the optical/UV continuum and the X- 
ray absorbing column as measured by XMM-Newton both appear to be dominated by the Galactic 
contribution, suggesting that the absorber associated with M87 has Nh < 10 20 cm" 2 . Recent 
Chandra X-ray measurements reported by Di Matteo et al. (2002) indicate that any absorbing 
matter associated with the nucleus has Nh < 3.2 x 10 20 cm" 2 . A consistent picture thus emerges 
in which the absorber has Nh ~ 10 19 — 10 20 cm" 2 . The only potential disagreement arises from the 
lower limit of Nh > 10 20 cm" 2 obtained above from the Mn II absorption lines. The Mn II A2577 
line is potentially contaminated by Galactic absorption in Fe II A2587, and Mn II A2594 suffers a 
similar problem with Galactic Mn II A2606, leaving the Mn II A2606 in M87 as the best option for 
measuring this ion; the measurement uncertainty for this feature is relatively large, however (la 
error of 26% ; see also Fig. 5c/), suggesting that it does not represent a significant inconsistency. 

A fundamental question concerning the absorber in M87 is its relation, if any, to the accretion- 
powered activity in the nucleus. LINERs commonly exhibit resonance-line absorption tracing mod- 
est column densities, and these absorbers characteristically display a low-ionization state that mir- 
rors the low ionization of the emission plasma in these sources (Shields et al. 2002, and references 
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therein). These absorbers are not particularly exotic in their properties, and a natural question is 
whether they have any connection to the active nucleus, or simply represent normal components of 
the interstellar medium of the host galaxy that happen to fall along our line of sight to the central 
UV source. The latter possibility is underscored by the similarity in detected ions and equivalent 
widths for the Galactic and M87 absorbers (Tables 4 and 5). 

While the nature of absorbers in LINERs in general remains ambiguous, in the case of M87 
there are indications that the absorbing material is, in fact, associated with the central accretion 
phenomenon. The offset in velocity between the absorber and the underlying galaxy, indicative 
of outflow of the absorbing matter, is highly reminiscent of absorption commonly associated with 
more luminous AGNs, which is often blueshifted with typical velocities of a few xlOO km s _1 (e.g, 
Crenshaw et al. 1999). In addition, Sankrit et al. (1999) concluded that the substantial strength 
of the Lya emission feature requires that the absorber cover the continuum source, but little if any 
of the emission-line clouds, if Nh > 10 18 cm" 2 , as seems to be the case. The most straightforward 
way of interpreting this finding geometrically is that the absorber has a distance from the central 
source and a size that are smaller than the typical dimensions of the emission clouds. Since the 
HST aperture employed by Sankrit et al. had a diameter of CK'26 ~ 20 pc, the size and scale of the 
absorber are probably much smaller than this dimension, which would be consistent with an origin 
in the central accretion structure or its environs. 

The absorption characteristics of M87 and possibly other LINERs may lead to important 
insights into the accretion process in these objects. The low ionization describing both absorp- 
tion and emission sets these sources apart from Seyfert nuclei and QSOs, and likely arises from 
a characteristically different mode of accretion at low accretion rates. In particular, physical sce- 
narios employing advection-dominated accretion flows ( ADAFs) , adiabatic inflow-outflow solutions 
(ADIOS), or other low radiative efficiency accretion structures are showing promise for interpreting 
these sources (see Quataert 2001 for a recent review). Substantial outflows of hot gas are predicted 
in some versions of these models (e.g., Blandford & Begelman 1999; Beckert 2000), and potentially 
some cooler material is accelerated in these winds, giving rise to the observed absorption. One 
quantity of interest for comparison with theoretical predictions is the rate of mass ejection M rep- 
resented by the M87 absorber. From simple geometrical considerations, M ~ fim p j "TV 'uA^rv , where 
[i is the mean atomic mass per H atom, m p is the proton mass, / is the global covering factor, r 
is the radial location of the absorber, and v is the outflow velocity. Assuming solar abundances 
(fi = 1.41; Dappen 2000), this can be expressed as 

M & 0.002/ ( 9 f H ) ( —) ( ^Meyr- 1 . (1) 

J Vl0 20 cm-V \lpc) Vl26kms-V ° V ' 

This result indicates that M is much less than the Bondi accretion rate for the central object, which 
is ~ 0.1 M© yr _1 (Di Matteo et al. 2002), but consistent with the idea that the absorbed gas may 
be material entrained in an outflow from the central accretion structure. Observations with higher 
spectral resolution will be necessary to establish the detailed kinematic and spatial structure of the 
absorbing gas in M87 and other LINERs. 
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4. Conclusions 

The present study of UV/optical emission and absorption in the nearby LINER M87 provides 
several new insights into the nature of this source. Within a radius of ~ 20 pc, emission-line 
ratios for the nucleus are well-described by photoionization of a multi-phase medium with a range 
of densities; shock models are less successful. The data also show a low-ionization absorber that 
evidently resides on small scales and is outflowing. These attributes, along with the famous radio 
jet, strongly resemble those of more luminous AGNs, and clearly point to accretion as the underlying 
power source responsible for the LINER behavior. 

M87 is, however, a good case for a composite system in terms of the ionization process generat- 
ing the optical LINER. The emission-line characteristics of the resolved nebula change significantly 
away from the center; the emission- line ratios (D97) as well as energetics arguments suggest that 
mechanical heating becomes important for this gas on scales of tens of parsecs. If this pattern 
is typical, larger apertures common in ground-based measurements will sample emission plasma 
that may trace a complex mix of energetic processes. Efforts to identify a single, or even domi- 
nant, ionization mechanism in these systems may thus be unproductive, although the question of 
whether alternatives to accretion power are important provides a good reason to conduct detailed 
investigations of other LINER galaxies that lack broad lines or other clear signatures of accretion. 
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Fig. 1. — STIS spectrum of the M87 nucleus. Undulations in the continuum redward of ~ 7700 A 
are artifacts of fringing in the CCD detector. 
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Fig. 2. — Optical and UV line ratio diagrams. Our measurements for the M87 nucleus are repre- 
sented by filled (STIS) and open (FOS) red circles, while green triangles represent measurements 
on the ionized disk (D97). Blue dotted lines show the loci of photoionization predictions, while 
magenta dashed lines indicate the ratios from the Dopita & Sutherland (1995) shock models. The 
low-density Am/7 sequence with parameters from Binette et al. (1996) is shown by the black dash- 
dot line with diamonds, while the high-density A m jj predictions described in §3.1 are represented 
by the black solid line with squares. (The square points for the models are suppressed in (b) for the 
sake of clarity.) For both A M /j sequences, successive points represent increases of A M /j by a factor 
of two; the smallest value is indicated. The asterisk indicates the composite model predictions 
for high-density A M /j = 0.004+ low-density IB. The effect of extinction (Ay = 1.0 or 10 mag) is 
shown by the arrow. 
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Fig. 3. — Same as for Figure 2, but for purely UV lines. 
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Fig. 4. — M87 nucleus spectrum (red line) and our model continuum fit (black line), consisting of 
a broken power-law modified by reddening. 
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Fig. 5. — Absorption features in the FOS and STIS spectra of M87 (solid histogram) and our model 
fits (dotted lines). 
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Table 1: HST STIS and FOS Observations 



STIS 






Spectral Resolution 1 




Grating 


Detector 


Coverage (A) 


Point Source 


Extended 


Exp. (seconds) 


G140L 


FUV-MAMA 


1150-1720 


190 


2500 


2590 


G230L 


NUV-MAMA 


1590-3160 


600 


4000 


2296 


G430L 


CCD 


2900-5710 


380 


1900 


900 


G750L 


CCD 


5270-10270 


390 


1900 


1000 


FOS 












Grating 


Detector 


Coverage (A) 


Point Source 


Extended 


Exp. (seconds) 


G160L 


FOS/BL 


1140-2510 


1300 


1300 


2300 


G270H 


FOS/RD 


2220-3275 


230 


230 


2300 


G400H 


FOS/RD 


3235-4780 


230 


230 


2300 


G570H 


FOS/RD 


4560-6820 


230 


230 


2300 



1 FWHM resolution in km s 1 . 
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Table 2: Observed Emission-Line Fluxes (STIS & FOS) 



Emission Line 


olio r lux 


r US r lux 


Lya 


14. U 1 ±u.oy 


on o/?_l_1 C\£. 

oU.obztz.yo 


U IV A1549 


O.o0±0.0fa 


0.99±U.28 


TT„ TT \ 1 d A C\ 

He II A 1640 


n /in i n c\n 

0.40±0.0o 


0.51±0.27 


ttt! \ i nnn 

Kj iiij A1909 


1.0o±U.13 


n no_i_n 07 
0.92±U.Z I 


P TTl 

^ 11J AZoZO 


l.UOznu.Uo 




Mg II AA2796, 2804 


1.12±0.09 


1.24±0.17 


[0 II] A3727 


2.47±0.11 


2.45±0.18 


[Ne III] A3869 


0.36±0.04 


0.79±0.09 


[S II] A4072 3 


1.52±0.08 


2.22±0.19 


H 7 4 


0.49±0.04 


0.52±0.08 


H/3 


1.00±0.05 


1.00±0.10 


[0 III] A4959 


0.61±0.03 


0.76±0.06 


[0 III] A5007 


1.84±0.08 


2.29±0.18 


[0 I] A6300 


1.62±0.08 


2.28±0.19 5 


[0 I] A6364 


0.53±0.03 




Ha, [N II] A6548,83 


12.33±0.49 


13.29±0.95 


[S II] A6717,6731 


3.10±0.13 


2.32±0.21 


[Ar III] A7165 6 


0.10±0.04 




[0 II] A7325 


0.87±0.08 





1 Normalized to H/3 = 6.45 x 10 

2 Normalized to H/3 = 4.03 x 10 

3 Blended with R5. 

4 Blended with [0 III] A4363. 

5 Includes [0 I] A6364. 

6 Uncertain. 



erg s cm , uncorrected for reddening, 
erg s _1 cm~ 2 , uncorrected for reddening. 
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Table 3: Photoionization Predictions 1 



Line 


high-n MB hig 


;h-n IB A M /i = 0.004 low-n IB Total 2 


Observed 3 


C IV A1549 


26.76 


0.05 


0.76 




0.00 


0.53 


0.60 


He II A 1640 


7.70 


0.08 


0.28 




0.70 


0.41 


0.40 


C III] A1909 


1.88 


1.26 


1.28 




0.10 


0.93 


1.03 


C II] A2326 


0.00 


1.11 


1.08 




0.66 


0.95 


1.05 


[0 II] A3727 


0.00 


0.05 


0.05 




6.48 


1.98 


9 47 

^ .4: 1 


[Ne III] A3869 


0.41 


1.82 


1.78 




0.96 


1.53 


0.36 


[0 III] A5007 


3.11 


1.91 


1.95 




0.32 


1.46 


1.84 


[0 I] A6300 


0.00 


2.33 


2.27 




1.47 


2.03 


1.62 


Ha + [N II] AA6548, 6584 


0.00 


3.80 


3.78 




9.81 


5.59 


12.33 


[S II] A6717,6731 


0.00 


0.26 


0.25 




6.91 


2.25 


2.32 


All fluxes expressed relative to H/3. 














Fluxes predicted for the composite, 3-component photoionization model. 




STIS measurements, uncorrected for reddening. 














Table 4: Absorption Lines (STIS, G140L) 








Line ID 


Kbs (A) 


EW(A) 


log(N cm 


2 ) 


Origin 




Si I A1255 1 


1259.31±0.16 


0.75±0.10 




14.35 




M87 




Si II A1260 2 


1264.49±0.33 


0.40±0.03 




12.42 




M87 




S I A1270 3 


1271.40±0.29 


0.33±0.12 




15.86 




Gal.(?) 




S I A1401 3 


1399.44±0.16 


0.79±0.10 




15.59 




Gal. 




Ni II A1413 3 


1417.96±0.38 


0.38±0.10 




15.55 




M87 




S I A1425 3 


1424.64±0.36 


0.47±0.14 




14.38 




Gal. 






1432.49±0.15 


0.39±0.09 






Unknown 






1435.13±0.41 


0.55±0.15 






Unknown 






1446.64±0.65 


0.59±0.16 






Unknown 




Ni II A1450 4 


1455.69±0.36 


0.50±0.16 




15.93 




M87 




Ni II A1468 5 


1468.55±0.79 


0.73±0.18 




16.30 




Gal.(?) 




Si II A1527 


1531.62±0.26 


0.70±0.19 




14.47 




M87 





1 Contributions from Galactic S I A1259.52 and Fe II A1260.53. 

2 Contributions from M87 S I A1259.52 and Fe II A1260.53. 

3 Multiple components. 

4 Contributions from M87 Co II A1448.01 and Galactic Ni II A1454.84. 

5 Contributions from Galactic Co II A1466.20 and Ni II A1467.26. 
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Table 5: Absorption Lines (FOS, G270H, G400H, G570H) 



Line ID 


Kbs (A) 




EW(A) 


log(N 


cm 2 ) 


Origin 


Fe II A2344 


2343. 


.23±0. 


.15 





.63±0. 


,11 


14 


.11 


Gal. 


Fe II A2344 


2352. 


.87±0. 


.39 





.59±0. 


.16 


14 


.08 


M87 


Fe II A2374 


2373. 


.86±0. 


.28 





.80±0. 


,24 


14 


.80 


Gal. 


Fe II A2383 


2382. 


.22±0. 


.22 


1 


.13±0. 


,15 


13 


.94 


Gal. 1 


Fe II A2383 


2391 


.09±0. 


.22 


1 


.33±0. 


,15 


13 


.04 


M87 


Mn II A2577 


2586. 


.21±0. 


.52 





.70±0. 


,48 


13 


.57 


M87 


Fe II A2587 


2586. 


.11±0. 


,39 





.39±0. 


35 


14 


.03 


Gal. 


Fe II A2587 


2595. 


.62±0. 


36 


1 


.21±0. 


,21 


14 


.59 


M87 2 


Fe II A2600 


2599 


.74±0. 


.24 





.73±0. 


.17 


13 


,78 


Gal. 


Fe II A2600 


2609 


.70±0. 


.26 


1 


.66±0. 


.16 


14 


.20 


M87 


Mn II A2606 


2606 


.03±0. 


.57 





.87±0. 


.21 


13 


.92 


Gal. 


Mn II A2606 


2617 


.19±0. 


.78 





.72±0. 


,19 


13 


.83 


M87 


Mg II A2796 


2795 


.91±0. 


.10 


1 


.28±0. 


.12 


13 


.58 


Gal. 3 


Mg II A2796 


2803 


.08±0. 


.11 


1 


.59±0. 


,40 


13 


.70 


M87 4 


Mg II A2803 


2813 


.61±0. 


.24 


1 


.91±0. 


,55 


14 


,05 


M87 


Mg I A2853 


2863 


.19±0. 


35 


1 


.27±0. 


,18 


12 


.08 


M87 


Ca II A3935 


3949 


.40±0. 


.26 


2 


.10±0. 


.20 


13 


.45 


M87 


Ca II A3970 


3984. 


.29±0. 


32 


1 


.83±0. 


,12 


13 


.80 


M87 


Na I A5891 


5909. 


.44±2. 


.15 





.61±0. 


,69 


12 


.48 


M87 


Na I A5897 


5915 


.01±1 


.07 


2 


.96±0. 


,72 


12 


.53 


M87 



1 Contaminated by M87 Fe II A 2376. 

2 Contaminated by Galactic Mn II A2594. 

3 Contaminated by Galactic Mn I A2795. 

4 Contaminated by Galactic Mg II A2803 and M87 Mn I A2795. 



